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ABSTRACT 

The origin of a double main-sequence (MS) in oj Centauri is explored. We have shown from theoret- 
ical calculations on the stellar evolution that the colors of MS stars are shifted to those of the observed 
blue MS if the surface layers are polluted by He-rich materials with the mass of ~ 0.1 Mq. Stars 
are supposed to be polluted through numerous encounters with the ejecta descended from massive 
asymptotic giant-branch (AGB) stars. Two populations of stars with different kinematics exception- 
ally observed in oj Gen indicate that kinematically cooler stars are more polluted through encounters 
with AGB ejecta than kinematically hotter ones because the accretion rate is inversely proportional 
to the cube of the relative velocity. We propose that both of these factors split the MS in w Gen. 
This theoretical scheme explains why only w Gen exhibit a double MS and matches the amount of 
He necessary to produce the blue MS with that supplied from massive AGB stars. Furthermore, we 
predict that even if globular clusters (GCs) possess only one generation of stars, the velocity dispersion 
of stars broaden the MS in the color-magnitude diagram as long as the GCs are massive enough to 
keep the AGB ejecta after the burst of star formation. This view explains the broad MS recently 
found in the GG NGG 2808 which exhibits no scatter in [Fe/H] and thus is likely to consist of a single 
generation of stars unlike the case of uj Gen. 

Subject headings: globular clusters: general — globular clusters: individual {uj Centauri, NGG 2808) 
— stars: abundances — stars: AGB and post- AGB — stars: evolution 



1. INTRODUCTION 

The color-magnitude diagram (CMD) obtained by re- 
cent Hubble Space Telescope and VLT observations have 
clearly revealed the existence of a double main-sequence 
(MS) in the Galac tic globular cluster (G G), lo Centauri 
(|Bedin et al.ll2004l: [Sollima et al.ll2006af) . In addition to 
the dominant red MS (rMS) population, the blue MS 
(bMS) occupies - 20 - 35% of the MS stars. Sur- 
prisingly, spectroscopic observations with GIRAFFE at 
the VLT have revealed that the rMS has a metallicity 
of [Fe/H]~ —1.6 with some scatter , whereas the bM S 

is more metal-rich: [Fe/H] 1.3 (jPiotto et al.ll2005f ). 

Thus, their metallicities do not explain this color dif- 
ference. A high content of He such as F ~ 0.35 — 0.4 
in the bMS stars is likely to be the only possible origin 
(|Bedin et al.|[200l : iNorris 2004). 

Red giants and horizontal-branch (HB) stars in 
OJ Ge n exhibit a large spread in their metal l icities 
(e.g.. INorris fc Da Costal [l99l ISuntzeff k Krai [l996l : 
ISollinia et al.l l2006b[ ). Especially, the metallicities of 
red giants spread over — 2 [Fe/H] ^ — 0.5, and ap- 
parently divide the red giant branch into three distinct 
sub-populations: the dominant (the fraction / ^ 80%) 
metal-poor ( [Fe/H] ;$ —1.4) population; intermediate- 
metallicity ( -1.4 < [Fe/H] < -- 0.8) less populous (/ - 
15%) one (|Norris et al.l |1996[ ): extremely metal-rich 
([Fe/H]~ —0.6) anomalous one, so called, RGB-a with 
/ ~ 5% of the en tire population (jPancino et al.l I2OOOI : 
iFerraro et al.|[2003 ). Furthermore, the difference in stel- 
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lar metallicities is linked to their kinematics: more metal- 
rich s tars are kinematic ally cooler than metal-poorer 
stars ()Norris et al.l 119971 ). Such complex features indi- 
cate multiple episodes of star formation in w Gen, which 
is exceptional among GCs because the other Galactic 
GCs exhibit essentially no metallicity dispersion proba- 
bly indicating a single generation of stars. The duration 
of star formation in uj Gen is inferred to be a few Gyrs 
from the observed elemental abund ance ratios of red gi- 
ants (|Tsuiimoto fc Shigevamal[2003[ ). In such a scheme of 
self-enrichment, the bMS stars should have been formed 
from the gas enriched by the former generations of stars 
evolved from the rMS. This picture naturally provides 
the condition that the amount (Ay ~ 0.15) of the extra 
He in the bMS stars needs to be less than the amount 
of He suppli ed from massive asy mptotic giant-branch 
AGB) stars (iD'Antona et al.l l2005D and/or massive stars 
Maeder fc Mevnetll2006[ ) evolved from the rMS popula- 
tion. If the bMS stars are entirely composed of He-rich 
gas, this co ndition has been sho^y n to be hardly satisfied 
in UJ Gen (jBekki fc NorrisI |2006[ ) and will be reviewed 
in t he following sec t ion. To avoid the short supply of 
He, 'Bekki fc N orrisI (|200 6^ proposed an alternative sce- 
nario in which He is supplied from an external sources in 
a dwarf galaxy whose nucleus became uj Gen afterward. 
This scenario also assumes that the bMS stars are formed 
from He-rich gases. 

In this letter, we propose a new scenario in which kine- 
matically cooler MS stars enhance their He contents in 
the surface layers by accreting He enriched matter when 
they encounter with AGB ejecta. Based on this scenario, 
we will shed light on the origin of the double MS as well 
as its uniqueness. As far as we know, uj Gen is the only 
object exhibiting a split in the MS. Similarly, NGG 2808 
is found to exhibit a broad MS in the CMD, a blue part 
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of which is claimed to be reproduced by a super He- rich 
popu l ation similar to the case of ui Cen (|D'Antona et al.l 
[20051: iLee et"all 120051) . However, in contrast to lo Cen, 
stars of NGC 2808 do not show a metallicity spread 
l|Carretta et al. ' '20041 Thus, the double MS should be 
caused by characteristics unique to uj Cen, while the ex- 
istence of the blue-shifted MS stars should be explained 
by common features shared with these two GCs. 

2. THE AMOUNT OF HE NECESSARY FOR BMS STARS 

Massive AGB stars that have experienced the second 
dredge-up are raised as a plausible production site of 
matter with F - 0.3 - 0.4 (|D'Antona et al.ll2005f) . Sup- 
pose that stars are born exclusively from the ejecta of 
these AGB stars, then the helium content of the stars 
will be close to 1^ ~ 0.4 and they will populate the bMS. 
On this occasion, we will check if the total mass of He 
in AGB ejecta is sufficient for making the observed num- 
ber of bMS stars in lo Cen. If we assume that the mass 
range of stars yielding Y ~ 0.3 — 0.4 is 3.5 — 8 Mq, 
we obtain the mass fraction /agb=0.08 of AGB stars 
yielding a high Y in the total stellar mass, adopting the 
Salpeter initial mass function (IMF). In terms of sup- 
ply and demand of He, the He mass supplied by such 
AGB stars is written as F x x /^ms x /agb x /ojocta 
where /ojocta denotes the mass fraction of AGB ejecta 
out of the progenitor star (=0.84). On the other hand, 
the He mass contained in the bMS stars could be writ- 
ten as F X X /bMS- Thus, these values combined 
with the observed mass fractions /bMS and /rMS of the 
bMS and rMS stars (=0.25, 0.75) suggest that He is 
in short supply in uj Cen, as shown by the relation: 
r X X /bMS > Y y. y. /^ms x /agb x /ojocta, 
where M,^ denotes the stellar mass of lo Cen. It has 
been pointed out that a variation in the IMF c an hardly 
change this magnitude relation (Bck ki fc Norr is 2006). 

Therefore, it is implausible that the bMS stars are as- 
sumed to be born from pure AGB ejecta. In the first 
place, little dilution of AGB ejecta with the interstel- 
lar matter before stars are formed, which is required to 
yield a high-He-content in the bMS stars, results in their 
metallicities almost the same as those of the rMS. In 
other words, the observed abundances of He and Fe in 
the bMS stars can not be consistently reproduced from 
this viewpoint. Even if stars were formed from pure AGB 
ejecta, there would be no reason why this happened only 
in w Cen. Besides massive AGB stars, [M > 10 Mq) 
with certain rotat ions are another candida te for yielding 
a high Y matter (jMaeder fc Mevnetl l200l l . This candi- 
date suffers from a similar shortage of He supply. 

3. SURFACE POLLUTION SCENARIO 

UJ Cen is the only GC in the Galaxy that consists 
of multiple generations of stars. These generations are 
also kinematically heterogeneous orris eiTanfTool . In 
this chemically and kinematically heterogeneous environ- 
ment, AGB stars eject their stellar end products includ- 
ing He in the later evolutionary stage, a part of which 
take the form of planetary nebulae (PNe). The size of 
the AGB ejecta reaches ~ 0.2 pc which is a typical size of 
PNe. As a consequence, such He-rich gases will occupy 
most of the volume of proto-tj Cen because GCs includ- 
ing LO Cen are very compact objects whose half radii are 
a few pc. Thus, kinematically cool metal-rich stars can 



accrete He enriched gas through encounters with these 
AGB ejecta. Then the accreted matter is mixed inside 
the surface convective layers and enhances the He con- 
tents to make their apparent colors to bluer. In con- 
trast, metal-poor stars moving faster cannot accrete He 
enriched gas as much as the metal-rich stars do because 
the accretion r ate is invers ely proportional to the cube 
of the velocity (|Bondilll952l ). 

This occurs exceptionally in a massive GC such as lo 
Cen because its deep gravitational potential can retain 
AGB ejecta. In this scenario the abundance of He is 
enhanced only in the outer layer of the metal-rich and 
kinematically cool MS stars. Therefore, it significantly 
reduces the amount of He needed to shift stellar colors 
to bluer. In order to explain the observed split of the 
MS t hat reaches down to Vfjoe 22 — 23 (jBedin et al.l 
l2004f ). the necessary mass Mace of matter with Y = 0.4 
accreted to each metal-rich star is found to be ~ 0.1 
Mq as shown in the next section. This is larger than the 
typical mass of convective envelope in low mass MS stars. 
Then the amount of ejecta from massive AGB stars of all 
generations of stars is found to be /agb x /ojocta x M^^ 
O.OTM^, which is greater than the mass supposed to be 
accreted by the bMS stars /bMS x Mace x M^/Mstar ~ 
0.03Mj^. Here Mgtar denotes the average mass of a star 
(=0.8 Mq), Mace the mass of the AGB ejecta accreted 
to the star (=0.1 Mq). 

There exists supporting evidence for the surface pol- 
lution scenario. RR Lyrae stars with their metallici- 
ties of [Fe/H]~ —1.2 corresponding to that of the bMS 
stars are found to ha ve no sign of the He enhancement 
(|Sollima et al.ll2006bD . It is a natural consequence from 
the surface pollution scenario because stars increase the 
masses of their surface convective layers and thus tend to 
erase the enhancement of He as they evolve from the MS. 
On the contrary, the observed fact that extreme (very 
hot) HB stars occupy a fraction of ~ 30% (jP'Cruz et al.l 
[2000 ) similar to bMS stars in to Cen may imply the link 
between bMS stars and extreme HB stars. This fe ature is 
expla ined in the fram ework of fully He-rich stars (^ Norrisl 
2004; ' Lee et al.l [20051 ). On the other hand, it has been 
suggested that the deep mixing in red giants is a promis- 
ing mechanism for abundance anomalies of red giants in 
GCs (iKraft I I1994D and hence it can be related to the ex- 
treme HB (;e.g.. lSweiga7nil997t[Suda fc Fuiimo"to1l2006D . 
Unravelling the second parameter problem is crucial to 
identify which mechanism is responsible for extreme HB 
stars in lo Cen. 

4. ACCRETION OF AGB EJECTA 

In this section, we will estimate the mass accreted 
by a kinematically cool star in lo Cen when it passes 
through He-rich AGB ejecta. A part of the envelopes 
of intermediate-mass stars are ionized during the mass- 
losing process at the post- AGB stage and are visible as 
PNe. Though the total mass ejected from an AGB star 
should be a few or several solar masses, observations of 
Galactic PNe have revealed only 0.2 Mq as a luminous 
mass on average. That is a so-called missing mass prob- 
lem in PNe. It has been suggested that most of the miss- 
ing mass forms molecular clumps. Such a dense form of 
AGB ejecta is possibly inclined to survive without dis- 
persing into the ISM for a prolonged time. Thus, we 
assume that most of AGB ejecta take a form of molecu- 
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lar clumps distributed over a region with the size rojccta 
corresponding to that of PNe and survive for a period 
of star formation in w C en that sho uld be ~ 2 — 4 Gyrs 
(|Tsuiimoto fc Shigevam a 2003; Stanf ord et aI1l2006f) . At 
the end of the star formation, the remaining gas and 
AGB ejecta are supposed to b e stripped during a pas- 
sage through the Galactic disk (jTsuiimoto fc Shigevamal 
|2003() . This scheme suggests that the surface metallic- 
ity of the bMS stars is conve rged to the averag e value of 
u Cen (< [Fe/H] >~ -1.4. iNorris et al.l [19961 ) because 
massive AGB cjccta descended from all populations can 
pollute the surfaces of stars. That is compatible with 
no detection of any metallicity spread in the bMS stars 
(jPiotto et al.| [2005"). 
The total mass Mace accreted to a star is expressed as 

Mace = dM^cc/dt ■ T^/t 

cnc ' ■^^cjecta 

where tone denotes the timescale for the encounter be- 
tween a star and AGB ejecta and v is their relative ve- 
locity. Each factor in the right hand side is the mass 
accretion rate, the total number of encounters per star, 
and the time for a star to traverse the region of AGB 
ejecta. First, the accretion rate due to the Bondi accre- 
tion is written as 

dMi^cc/dt ^ nrgPcjcctaV, (2) 

where rg{= GMstar/v^) denotes the accretion radius of a 
star and Pojoeta is the density of the AGB ejecta averaged 
over the region with a radius of rejeeta- To evaluate the 
accreted mass, we adopt input values of Afstar = 0.8 Mq, 
V = 10 km s~^, and rejeeta = 0.2 pc together with the 
mass of AGB ejecta — 4.2 Mq which is deduced from the 
average mass of 3.5-8 Mq AGB stars (= 5 Mq) and the 
mass of a white dwarf (= 0.8 Mq). Then, we obtain 

dM^^Jdt - 3.3 X lO^-* g/s. (3) 

Next, the timescale <enc is obtained by tone = 
(»^cjeetaO-ejoetaw)~\ whcrc Hejeeta is the number density 
of AGB stars that can supply He in the proto-w Cen 
whose radius is , and (Tejeeta is the cross section of each 
ejecta given by T^r'^jccta- Here we assume the constant 
'T-cjeeta for all period of = 2 Gyr. The total num- 
ber of AGB stars is estimated to be 3.3 x lO** in the 
proto-tj Cen with the mass of 2 x 10^ Mq, assuming the 
Salpeter IMF and the mass range of 3.5-8 Mq relevant 
to the major He contributor. This number together with 
r^j = 6 pc which corresponds to the present half radius 
of (jj Cen yields riejeeta 1.4 x 10~^'*cm~'^. Then, we 
obtain tenc ~ 2.2 x 10'^ yr, which gives 9.2 x 10^ encoun- 
ters in total for the period of t^. Finally, combined with 

2rejoeta/'y ~ 1-2 X lO^^g^ dcduCC Maee ^ 3.6 X 10^^ g - 

0.18 Mq. The derived accretion mass is comparable to 
the mass of the convective envelope of a star with the 
mass of ^^0.4 Mq. This amount of the accreted mass is 
sufficient to realize a double MS which continues from 
the MS turn-off down to Vqog -- 22 - 23 in the CMD. 

5. COLORS OF STARS POLLUTED BY HE-RICH MATTER 

5.1. Isochrone 

We calculate the stellar evolution of low-mass stars 
with the masses of 0.3 to 0.85 Mq, the surfaces of which 
are polluted by a y = 0.4 matter. The mass accreted to a 
star is assumed to be 0.1 Mq. For simplicity, we neglect 
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Fig. 1. — Theoretical CMDs of the main sequence of w Cen cal- 
culated for the assembly of stars with y = 0.24 (solid curve) and 
of those stars that accrete the Y = 0.4 matter by the amount of 
O.IM0 (dashed curve). For reference, the observed points (filled 
triangles) and ranges (crosses and bars) are shown at four vi- 
sual magnitudes of 20, 21, 22, and 23 mag for bMS and rMS, 
respectively. Thes e observed data are extracted from Fig. Ic of 
IBedin et all ll200i) . 



a duration of mass accretion and thus place the y=0.4 
matter on the surface layer of a star at the start of calcu- 
lations. These results are compared with those for stars 
without a pollution. The stellar evo lution p rogram and 
input physics are the same as in Suda et all (12004) . The 
initial abundances are set to be ¥=0.24 and [Fe/H] = 
— 1.6, —1.3 for two populations, and y=0.40 and [Fe/H] 
= — 1.3 for the accreted matter. All models are assumed 
to be a-enhanced so that [a/Fe]=0.4. These computa- 
tions do not include any special treatment for convection. 
The border of convective zones are determined by the 
Schwarzschild criterion with mixing length parameter of 
1.5. Isochrones are obtained through the conversion of 
calculated effective temperatures and luminosities into 
the color Vgoe ~ and the magnitude Vgoe of the HST 
filters according to the recipe of Hol tzman et al.l (|1995f) . 
Extinctions are also c alculated from the values given in 
iHoltzman et al.l ()1995l ). Here we adopt the distance mod- 
ulus (m - M)=13.4 and the reddening E{B - V)=0.15. 
The ages of the two populations are assumed to be 13 
Gyr with A(Age)=0 since red giants corresponding to 
both populations exhibit Type II supernova-like abun- 
dance patterns suggesting that the age difference should 
be less than ~ 1 Gyr and could be ignored. 

Figure 1 shows the Vgoe — -^814 vs. Vgoe diagram pre- 
dicted for F=0.24 stars (solid curve) and stars with a 
O.IM0 accretion (dashed curve), together with yard- 
sticks and filled triangles indicating the ob served loca- 
tions of the rMS and the bMS, respectively (|Bedin et al.l 
l2004l ) . The isochrone of stars with the accretion lies close 
to that of the observed bMS. In our models, the MS split 
diminishes at around Vgoe ~ 23, which is compatible with 
the observation. 
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Fig. 2. — The distribution functions with respect to the color 
^606 ~ ^814 of main sequence stars accreting He-rich AGB ejecta. 
The dashed line shows the distribution function of main sequence 
stars with [Fe/H]= —1.6 and the velocity dispersion of err = 60 
km/s and the dash-dotted line the distribution function of main 
sequence stars with [Fe/H]= —1.3 and cr^, = 10 km/s. The sum 
of these two populations are shown by the solid line. The number 
fraction of the blue main sequence stars is assumed to be 35%. The 
total number of stars are arbitrary. 



5.2. Double main sequence 

The actual color shift of stars due to the surface pollu- 
tion discussed in the preceding sections should be varied 
because each star moving at a different speed accretes 
He-rich gas with a different mass accretion rate. We 
have derived the distribution function of such stars with 
respect to the color (Vgoe ~ -^814) using the formula 



dN 



d (Veoe ^ -^814) 



dN dv 



dAL. 



dv dAfacc d {Veoe - hu) 



(4) 



where dN/dv denotes the distribution function of stars 
with respect to the velocities. Here we assume the 
isothermal distributions with constant velocity disper- 
sions CTr for the rMS stars and at for the bMS stars. The 
accretion rate due to the Bondi accretion specifies the 
functional form of dw/dAfacc , which is proportional to the 
biquadratic of the velocity. The third factor is derived 
from the stellar evolution calculations (see §5.1|) for stars 
with the same total mass but with different masses of the 
accreted He-rich matter. Since the derived distribution 
is for stars with the same total mass, this will give the 
color distribution of stars with the approximately same 
Veo6i= 19.3 ±0.2). 

Though the colors of the rMS stars are blue-shifted as 
well as the bMS stars, the amount of the shift is signif- 
icantly larger for the bMS. This is because the velocity 
dispersion db for stars with [Fe/H]~ —1.3 is smaller than 
(Tr for metal-poor stars. A larger velocity dispersion re- 
duces the mass of accreted gas and hence the amount of 
the color shift due to the He enhancement. For exam- 
ple, the distribution with respect to the color for stars 
with the mass of 0.77 Mq is presented in Figure ??. Here 
the velocity dispersions are assumed to be CTi- = 60 km 
and (Tb = 10 km s^^. The distribution has a clear 



gap between the rMS and bMS stars. Although the 
value of ffr seems extremely large in comparison with 
the present velocity dispersion in u; Cen, an elemental 
abundance a nalysis is suggestive of ~ 50 km/s in the 
proto-tj Cen (Tsuii moto fc Shig evama ''2003V The distri- 
bution function of the rMS stars presented in Figure ?? 
has a too sharp peak as compared with the observations 
([Sollima et al.ll2006al ) because our calculations ignore the 
scatter in metallicity observed for the rMS stars. This 
scatter will broaden the rMS to redder parts of the CMD. 

Finally, it should be of note that our scenario implies 
that some GCs with a single generation of stars hold a 
broad MS rather than a double MS as long as these GCs 
are massive enough to retain AGB ejecta in their deep 
gravitational potentials and stars can pass through AGB 
ejecta many times. 



6. CONCLUSIONS AND DISCUSSION 

We propose a new mechanism to produce a double MS 
in the CMD of uj Cen, focusing attention on the fact 
that high He abundances in the stellar surface layers in- 
crease the surface temperatures. On the other hand, all 
the other schemes proposed so far assume that stars in 
the bMS have high He abundances not only in the surface 
layers but in the entire volumes. However, these schemes 
confront severe difficulties; (i) the amount of He supplied 
from AGB stars and/or supernovae evolved from the red 
MS population would not be sufficient to produce all the 
bMS stars, (ii) the observed abundances of He and Fe in 
the bMS stars can not be consistently reproduced, and 
(iii) there is no reason why this rare phenomenon oc- 
curred only in u! Cen. These critical issues are settled 
in our scenario as follows. Kinematically cool stars ac- 
crete a sufficient amount of AGB ejecta to become the 
bMS stars. The surface abundance of Fe of these stars 
converges to the average value of u Cen. In contrast, 
kinematically hot stars accrete too little AGB ejecta to 
significantly change the colors. Detailed elemental in- 
formation besides He or Fe for individual MS stars in uj 
Cen will give a deeper insight into the origin of blue MS 
stars. In fact, a preliminary result on the average C and 
N abundances of MS star s seems to pose a problematic 
challenge to our scenario (jPiotto et al.|[200l . 

Generally, a GC undergoes a single burst of star forma- 
tion because subsequent supernova explosions heat and 
evaporate the interstellar gas from the GC and thereby 
stars never accrete He-rich gases. On the other hand, 
some GCs are expected to be massive enough to re- 
tain the interstellar gas (|Tsuiimoto fc Shigevamall2003f) 
as well as the AGB ejecta and allow stars to accrete 
He-rich gas through encounters with the AGB ejecta. 
Among these GCs, w Cen is an exceptional case in the 
sense that it has multiple stellar populations with dif- 
ferent kinematics. Here we have shown that this envi- 
ronment surely results in a split in MS. The mechanism 
presented here predicts that other massive GCs with suf- 
ficiently deep gravitational potentials to hold AGB ejecta 
host bluer MS stars that accreted He-rich gas even when 
no kinematically different population exists. Each star 
in a GC has a different velocity obeying a certain dis- 
tribution and thus accretes a different amount of AGB 
ejecta. NGC 2808 possibly belongs to these massive GCs 
without multiple kinematic populations. 
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